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ABSTRACT 

We analyse the effect of Hall currents on the Convective Heat 
Transfer flow of a viscous electrically conducting fluid in a 
vertical channel bounded by the flat plates at X = + L. Which are 
maintained at Non uniform Temperature in the presence of Heat 
sources taking the slope S of the boundary temperature the Non- 
linear coupled equations governing the flow and Heat Transfer are 
solved by employing a regular perturbation technique the velocity 
and temperature distributions are analysed for different values of 
the Hartmann Number M, Hallcurrent parameter m, Grashoff 
Number G, and Radiation parameter N[. The rate of Heat 
Transfer has been calculated numerically for different values of 
the governing parameters. 

Keywords: Heat Transfer, Hall effects, Heat sources, Vertical 
channel. 



1. INTRODUCTION 

The magneto hydrodynamic heat 
transfer has gained significance in recent 
times owing to its applications in recent 
advancement of space technology. The 
process of free convection as a mode of heat 
transfer has wide applications in the fields of 
chemical Engineering, Aeronautics and 
Nuclear power generation. Natural conve- 



ction in porous media with internal heat 
generation is of interest in such situations as 
post accident heat removal in nuclear power 
reactors and the geothermal problems arising 
during the storage of nuclear waste in the 
earth 2 . Rajeswara Rao 6 has analysed the 
combined free and forced convective flow of 
an electrically conducting, viscous, 
incompressible fluid confined in a vertical 
channel whose boundaries are maintained at 
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a non-uniform temperature He has not 
considered induced magnetic field into 
account Taking induced magnetic field into 
account Ram Chandra 7 has discussed the 
free and forced convective flow of an 
electrically conducting fluid in a vertical 
channel whose walls are maintained at non- 
uniform temperature Ravindra 8 has 
investigated the natural convective flow and 
heat transfer through a porous medium in a 
vertical channel maintained at Non-uniform 
temperature with constant heat sources. 

The unsteady flow of a rotating 
viscous fluid in a channel maintained by 
non-tensional oscillations of one or both the 
boundaries has been studied by several 
authors to analyse the growth and 
development of boundary layers associated 
with geothermal flows for possible 
applications in geophysical fluid 
dynamics 2 ' 3 ' 4 ' 5 ' 9 . Later Singh et al. m studied 
free convection in MHD flow of a rotating 
viscous liquid in porous media. Singh et al. 1 1 
have also studied free convective MHD flow 
of a rotating viscous fluid in a porous 
medium past an infinite vertical porous 
plate. However, in a partially ionized gas, 
there occurs a Hall current when the 
strength of the impressed magnetic field is 
very strong. 

2. FORMULATION 

We consider the steady flow of an 
incompressible, viscous, electrically conduc- 
ting fluid confined in a vertical channel 
bounded by two flat walls under the 
influence of an inclined magnetic field of 
intensity Ho lying in the plane (y-z).The 



magnetic field is inclined at an angle a to 
the axial direction k and hence its 

components are (0,H Q Sin(a),H Q Cos(a)). 

The walls are maintained at non-uniform 
temperature.In view of the non-uniform 
boundary temperature the velocity field has 
components(u,0,w)The magnetic field in the 
presence of fluid flow induces the current 
(J x ,0,J„). We choose a rectangular 

cartesian co-ordinate system 0(x,y,z) with z- 
axis in the vertical direction and the walls at 
x = ±L. 

When the strength of the magnetic 
field is very large we include the Hall 
current so that the generalized Ohm' s law is 
modified to 

J + 0) e T e JxH = <J(E + jU e qxH) (1) 

where q is the velocity vector.H is the 
magnetic field intensity vector.E is the 
electric fieldj is the current density vector, 
0) e is the cyclotron frequency, T e is the 
electron collision time, a is the fluid 
conductivity and ^is the magnetic 

permeability. 

Neglecting the electron pressure 
gradient,ion-slip and thermo-electric effects 
ans assuming the electric field E=0,equation 
(2.6) reduces 

J x - m H 0 J z Sin(a) = -oju e H 0 wSin(a) 

(2) 

J z + m H 0 J x Sin(a) = <Jju e H Q u Sin(a) 

(3) 

where m= 0) e T e is the Hall parameter. 
On solving equations (2) & (3) we obtain 
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J, = 



J, = 



0[i e H a Sin(a) 
\ K l+m 1 HlSin 1 (a) 



(J{i e H a Sin(a) 
1+m 2 HlSin\a) 



(mH 0 Sin(a)-w) (4) 



a = 



(Heat Source Parameter) 



where u,w are the velocity components 
along x and z directions respectively, 



(u + mH 0 wSin(a)) (5) The corresponding boundary conditions are 

W(f)-¥(-f)^ 



On introducing the following non- dlf/ 



dimensional variables 



, , (x,z) , W T-T e 



= <d,^ L = Q,0 = y(dz) at x = -L 
oz ox 



-^ = 0,-^ = 0,0 = ;K&) at x = +L 
dz ax 



The equations of momentum and energy in 

the presence of heat generating sources in 3. METHOD OF SOLUTION 
the non-dimensional form are 



vy-M! 2 vy+ 



RAdx. 



R 



PR 



dz )[ dx 



dx 



dz 



J J 

(6) 

mymmy^ m 



where 

f 

G = 



3 \ 



M 2 = 
R = 
P 



v 



(Grashof Number) 



V ' J 



Introduce the transformation such 

that z = Sz, t~ = S^- 
oz oz 

~\ ~\ 

Then— ~0(S)^ — -0(1) 
oz oz 

For small values of 8 «l,the flow develops 
slowly with axial gradient of order 8 and 

hence we take — -0(1). 
dz 

Using the above transformation the 
equations governing equations reduces to 



V 



(Hartman Number) 



F V -m, 2 fV + 



d 0 

3 x 



qL 



V v J 



V K f J 



(Reynolds Number) 
(Prandtl Number) 



8 R 



SPR 



d \/f 8 ( F V ) _ 9 y d ( F 2 y ) 
3 Z~~ d x d x d~ 



(8) 



F l e-ae 



(9) 
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F 2 =^ + S 2 1_ 
dx dz 



Assuming the slope 5 of the wavy 
boundary to be small we take 

y/{x,z) = W 0 ( X > y) + Syf^z) + 
S 2 y/ 2 {x,z) + 



0(x,z) = 0 o (x,z) + S0 1 (x,z) + 
S 2 d 2 {x,z) + 



(10) 



Substituting (10) in equations 
(8)&(9) and equating the like powers of 8 
the equations and the respective boundary 
conditions to the zeroth order are 



dy 



j- - (M | ) - 



dy 



R A dz 



d v 0 

dy 



3> 0 ] (dy/A( a> 0 



^ dz' J V 3z Jydxdy 



with 



(15) 



^(+l)-^(-l) = 0 

^=0,^=0, 0 1= O atx = -\ 
dy dz 

iZL =0 ,^=0, ffl =0 atx = + l 

dy dz 



Where 



N 7 = 3Nl , a,= aN 7 , P. =PN ? 
2 3N l +4 1 21 2 



(16) 



dy 2 



a x 0 Q = 0 



— L±-(M 2 ) — — 

^ dy dy J 



with 

j/ 0 (+l)-j/ 0 (-l) = l 
n 3 ^o 



: 0, 
n 



9y 



= 0, 



dz 



^d6p 

V dz , 



(11) 



(12) 



= 0, e a =y(z) at x = -\ (13) 
= 0, 0 a = y{z) at x = +1 



and to the first order are 

r d 2 0, 



ad, 



= PR 



[dy 

3x A 3z J I 3^ ){ dy 



(14) 



4. RATE OF HEAT TRANSFER 

The rate of Heat transfer (Nusselt 
Number) on the walls has been calculated 
using the formula 



Nu=- 



1 



(o m -e w ) 



dy 



m W ' V J J y=±\ 



where 0 m =0.5 j&dy 



and the corresponding expressions are 

(MO y= -i = 1 (-a 92 + Sa 94 ) 
yPm ~ 1) 

(Nu) y=+l =^-(a 92 + Sa 93 ) 



Journal of Pure Applied and Industrial Physics Vol.1, Issue 3, 30 April, 2011, Pages (162-211) 



N.B.V.R.D. Prasad, et al, J. Pure Appl. & Ind. Phys. Vol.1 (3), 170-177 (2011) 174 



5. ANALYSIS OF THE NUMERICAL 
RESULTS 

In this analysis we investigate the 
effect of Hall currents on Convective Heat 
Transfer flow of a Viscous Incompressible 
electrically conducting fluid in a Non- 
Uniformly vertical channel. We take the 
Prandtl Number P=0.71 and 8 = 0.01. The 
Velocity components u,w and The Non- 
dimensional Temperature distribution 9 are 
shown in Figs. 1-6. For different values of 
M,m, and Nr 

The Secondary velocity(u) which is 
due to the Non-uniform Temperature on the 
boundary is depicted in Figs.l and 2 for 
different parametric values. It is noticed 
from Fig. 1 represents the variation of u with 
M and m it is found that the secondary 
velocity retards with increase in Hall 
parameter m while it enhances with increase 
in the Hall parameter m in the entire flow 
region. 

Fig.2 shows that the variation of u 
with Radiation Parameter Ni shows that an 
increase in Ni results in an enhancement in 

H ■ 

The axial velocity w is shown in 
Figs. 3 and 4. The variation of w with 
Hartmann number M and Hall parameter m 
is exhibited in Fig.3. It is observed that 

higher the lorentz force smaller \w\ in the 

flow region also |w| enhances with increase 

in the Hall parameter m. The effect of 
radiation parameter Ni on the axial flow w is 



shown in FigAIt is found that an increase 

in Ni<5, \w\ depreciates in the left region 

and enhances in the right region and for 

higher N^IO. \w\ depreciates in the flow 

region except in the narrow region adjacent 
the lower boundary y=-l and this region 
reduces with increase in Ni. 

The Non-dimensional Temperature 
distribution (9) is exhibited in Figs. 5 and 6 
for different values of M,m, and Ni.lt is 
found that the temperature is positive for all 
parametric values variation. The variation of 
9 with M and m shows that higher the 
lorentz force larger the actual temperature 
and for further increase in the lorentz force 
smaller the actual temperature in the entire 
flow region.An increase in the Hall 
parameter m<1.5 depreciates the actual 
temperature and for further higher 
m>2.5.We notice an enhancement in the 
vicinity of the boundary y = ±l and 
depreciates marginally in the central region 
of the flow (Fig.5). The variation of 9 with 
Radiation parameter Ni .It is found that 9 
reduces with increase in Ni the inclusion of 
the radiative Heat Transfer leads to the 
remarkable depreciation in the Non- 
dimensional Temperature 9 (Fig. 6). 

The average Nusselt Number (Nu) 
which represents the which measures the 
rate of Heat Transfer at y= ± 1 is shown in 
Tables. 1-4 for different variations of 
G,M,m,a , and R.It is found that the rate of 
Heat Transfer depreciates with increase in 
the G>0 and enhances with G<0 at both the 
boundaries.The variation of Nu with 
Hartmann Number M shows that the higher 

the lorentz force larges \Nu\ at y=+l and 
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smaller \Nu\ at y=-l in the heating case 
while in the cooling case smaller |A/w| for 
M<15 and larger \Nu\ with M>20 at y=+l 
and at y=-l smaller \Nu\ an increase in the 
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Hall parameter m enhances |A/w| in the 

heating case and in the cooling case it 
reduces with m<1.5 and enhances with 
m>2.5.Also it enhances with increase in the 
Heat source parameter a and Reynolds 
Number R. (Tables. 1 and 2). 



Table . 1 Average Nusselt Number (Nu) at y=+l P=0.71, x=ji/4, (Xi=0.3, Ni=1.5 



G 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


10 3 


0.5514 


0.5822 


0.7438 


0.6635 


0.6915 


2.3829 


-20.9774 


0.6228 


0.7665 


2X10 3 


0.5097 


0.5593 


0.7272 


0.7312 


0.7180 


2.4844 


-19.5531 


0.6015 


0.7561 


-10 3 


0.6364 


0.6288 


0.7772 


0.5379 


0.6405 


2.2118 


13.4060 


0.6657 


0.7872 


-2X1 0 3 


0.6797 


0.6524 


0.7941 


0.4797 


0.6161 


2.1390 


8.7229 


0.6872 


0.7974 



Table . 2 Average Nusselt Number (Nu) at y=-l p=0.71, x= it /4, 0Ci=0.3, Ni=1.5 



G 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


10 3 


-0.4734 


-0.3979 


-0.2447 


-0.5883 


-0.6207 


-0.21886 


13.5318 


-0.5378 


-0.6635 


2X10 3 


-0.4323 


-0.3757 


-0.2295 


-0.6541 


-0.6463 


-2.2773 


18.3357 


-0.5164 


-0.6528 


-10 3 


-0.5574 


-0.4428 


-0.2753 


-0.4664 


-0.5718 


-2.0373 


-12.5995 


-0.5805 


-0.6848 


-2X10 3 


-0.6002 


-0.4656 


-0.2908 


-0.4099 


-0.5484 


-1.9732 


-8.2040 


-0.6019 


-0.6953 



M 


10 


15 


20 


10 


10 


10 


10 


10 


10 


m 


0.5 


0.5 


0.5 


1.5 


2.5 


0.5 


0.5 


0.5 


0.5 


a 


2 


2 


2 


2 


2 


4 


6 


2 


2 


R 


35 


35 


35 


35 


35 


35 


35 


70 


140 
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